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Stylophora pistillata is a widely used coral "lab-rat" species with highly variable morphology and a broad 
biogeographic range (Red Sea to western central Pacific). Here we show, by analysing Cytochorme Oxidase I 
sequences, from 241 samples across this range, that this taxon in fact comprises four deeply divergent clades 
corresponding to the Pacific- Western Australia, Chagos-Madagascar-South Africa, Gulf of 
Aden-Zanzibar-Madagascar, and Red Sea-Persian/ Arabian Gulf-Kenya. On the basis of the fossil record of 
Stylophora, these four clades diverged from one another 51.5-29.6 Mya, i.e., long before the closure of the 
Tethyan connection between the tropical Indo-West Pacific and Atlantic in the early Miocene (16-24 Mya) 
and should be recognised as four distinct species. These findings have implications for comparative 
ecological and/or physiological studies carried out using Stylophora pistillata as a model species, and 
highlight the fact that phenotypic plasticity, thought to be common in scleractinian corals, can mask 
significant genetic variation. 



DNA barcoding, usually based on the mitochondrial COI fragment 1 , has been extensively used to discrim- 
inate between closely related species, to identify new, cryptic or invasive species, and to assess biodiversity 
across many animal phyla 1 . The rule-of- thumb in DNA barcoding is that interspecific COI divergence is 
generally > 2%, whereas intraspecific variation is < 1 % 1 . However, based on 90 species from 44 genera belonging 
to 14 families, this criterion has been suggested as not being appropriate for most scleractinian corals 2 
(or anthozoans in general) because this region evolves very slowly in these organisms and consequently both 
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inter- and intra- specific variation are extremely low. We re-evaluated 
this conclusion by examining COI divergence and phylogeny of the 
"lab-rat" scleractinian coral, Stylophora pistillata, a species, which 
has been the focus of coral research over the last four decades. S. 
pistillata is a widely-distributed coral species in the Indo-Pacific, with 
numerous morphological variations (morphotypes) across different 
habitats, depths, and geographic regions 3 . Revisiting relationships 
between these morphotypes is necessary in order to ensure the cor- 
rect taxonomy for comparative studies based on this common coral 
species 4 . 

Results 

We obtained COI sequences from 241 S. pistillata colonies collected 
from 34 locations across most of its range (Western Pacific Ocean, 
Indian Ocean, Red Sea and Persian/ Arabian Gulf, Fig. 1). Interspe- 
cific variation in COI sequences was one of the highest documented 
among scleractinian genera (sequence diversity (p-distance) = 
0.01245 ± 0.00434), only surpassed by Pontes (0.02686 ± 0.0030) 
and Siderastrea (0.0141 ± 0.0027), which are globally distributed 
genera consisting of multiple morphologically divergent species 
(Fig. SI). High variation of the COI gene among these 3 genera sug- 
gests very early divergence of species within these genera 5 . Porites 
and Siderastrea are extant in the Atlantic, Pacific and Indian Oceans, 
whereas Stylophora is extinct in the Atlantic and only survived in the 
Indo-Pacific after the Plio-Pleistocene 5 . 



Phylogenetic construction based on COI sequences found 4 dis- 
tinct clades of S. pistillata (Fig. 2A). Clade 1 is comprised of speci- 
mens from the entire Pacific region and East Indian Ocean, ranging 
from New Caledonia in the southwest Pacific, Shfkoku (Japan) in the 
northwest Pacific, to western Australia in the Indian Ocean (Pacific- 
Western Australia clade, PWA). Clade 2 includes specimens from 
Chagos, St. Brandon's Island (Mauritius), La Reunion, Zanzibar, 
South Africa, and Madagascar (Chagos-Madagascar-South Africa 
clade, CMS A). Clade 3 includes specimens from the east African 
coast, the west coast of Madagascar and the Gulf of Aden (Gulf of 
Aden-Zanzibar-Madagascar clade, AGZM). Clade 4 includes the 
northwest Indian Ocean (Red Sea-Persian/ Arabian Gulf and 
Kenya, RSPAGK). The genus Seriatopora formed a fifth clade 
embedded at the base of clades 1 and 2. 

Intra-clade p-distance based on COI (Table SI) was considerably 
low (clade 1: 0.00017 ± 0.00009, clade 2: 0.00000 ± 0.00000, clade 3: 
0.00460 ± 0.00193, and clade 4: 0.00017 ± 0.00012). Higher intra- 
clade distance for clade 3 was due to the presence of specimens from 
Yemen that made a separate sub-clade into clade 3 (Fig. 2A). Inter- 
cede p-distance (Table SI) was always higher than intra-clade p- 
distance and ranged from 0.01274 ± 0.00436 (clade 3 vs clade 4) to 
0.01592 ± 0.00492 (clade 1 vs clade 3), with the exception for clade 1 
vs clade 2 (0.00378 ± 0.00248) and this could be due to the recent 
split between clade 1 and 2 compared to that between other clades 
(Fig. 3). The results form p-distance analyses also suggested the 
presence of 4 clades of Stylophora (Table SI). 




••• 

Clade 1 



Pacific Ocean 
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Figure 1 | Map showing the locations from which Stylophora pistillata samples were obtained (shown as red dots) . The division of four regions based on 
the four clades of Stylophora is marked in colours. Clade 1 - Pacific-Western Australia clade (Pink), Clade 2 - Chagos-Madagascar-South Africa clade 
(Green), Clade 3 - Gulf of Aden-Zanzibar-Madagascar clade (Blue) and Clade 4 - Red Sea-Persian/Arabian Gulf-Kenya clade (Brown). The samples were 
obtained from following locations; Japan (Kochi, Okinawa), Taiwan (Beitou, Penghu, Green Island, Renting), Taiping Island, Tioman, Indonesia, 
Australia (Lord Howe, Lizard Island, Western Australia), New Caledonia, The Chagos Archipelago, Re Union, Mauritius (St. Brandon's), Madagascar 
(North and South), South Africa, Renya (Ranamai), Zanzibar, Yemen, Dijibouti, Saudi Arabia (Arabian Gulf and Gulf of Aqaba), Oman, Iran and Red Sea 
(Eilat). The dotted green lines represent hypothetical extension of clade 4 and clade 2 to Pakistan-North India (Gulf of Rachhch) and Laccadives (west 
coast of India), south India and Sri Lanka. Area covered by blue line in the inset map shows the distribution of Stylophora pistillata (modified from the 
map in www.coralgeographic.com). The maps were drawn using the software Magic Maps Ver. 1.4.3 and Adobe Illustrator CS5 (Macintosh version). 
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Figure 2 | Phylogenetic analysis of 241 COI sequences of Stylophora pistillata. (A) Bootstrap values are based on Baysian/ML/NJ analyses. The numbers 
in the brackets stand for the total number of samples obtained form each location. Morphotypes and corallite structure of the four clades of Stylophora. 
(B) Based on SEM photos, in clade 1, six regular primary septa that are confluent to the center, fusing with the columella and corallites, are separated by 
the ccenosteum covered with sharp, fine and thin spinules. In clade 2, the arrangement of primary septa is similar to clade 1, but the spinules in the 
ccenosteum are more rounded, and fused corallites appear immersed in the ccenosteum. In clade 3, there are six poorly-developed primary septa with a 
rudimental second cycle of septa visible in some corallites. Septa remain separated and not confluent with the center. A columella is absent and the 
corallites are separated by poorly developed walls that raise to the surface of the ccenosteum. Spinules are rounded and poorly developed. Finally, in clade 
4, there is one cycle of six poorly-developed septa that is continuous with the hoods in the ccenosteum; the columella is well developed but is not fused with 
septa. Hoods are rounded and regularly distributed. Seriatopora morphotype and corallite structure is shown between clade 2 and clade 3. 



The four Stylophora clades were also supported by phylogenies 
based on nuclear ITS and calmodulin intron 3 (CAD3) (Fig. S2, Table 
S2). The deletion in CAD3 and duplication of tRNAW support the 
separation of the clades. In case of CAD3, complete deletion occurred 
in the samples from Pacific (clade 1), mixture of deletion and non- 
deletion in the central Indian Ocean (clade 2 and 3) and complete 
non-deletion in the Red Sea (clade 4). Similar results were also seen 
from the analysis of tRNAW duplication events (Table S2, Fig. S3). 



Comparison of tRNAW sequences from four clades showed duplica- 
tion event only in clade 1 and clade 2 (Table S2). Sequence results of 
tRNAW were also supported by restricted fragment length poly- 
morphism gel pattern using Rsal restriction enzyme (Fig. S3). 
While clade 1, 2 and 4 showed distinct patterns, clade 3 had mixture 
of clade 2 and 4. The results from the analysis of the Symbiodinium 
clade association showed that Stylophora in the pacific (clade 1) 
associated with Symbiodinium clade C and in the Red Sea (clade 4) 
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Figure 3 | Results from molecular dating on the phylogenetic tree of Stylophora. Five major ancestral nodes (A, B, C, D and E) to the clade relating 
to the four described clades of Stylophora were labeled to represent the major divergence events in the evolutionary history of Stylophora. Numbers on 
nodes are divergences time for each clades. The table shows the divergence time of each clade on the phylogenetic tree of Stylophora. Time of first 
appearance of Stylophora was used as reference time on the Bayesian evaluation of divergence on Beast (65-70 mya) . Results of molecular dating are listed 
in table which includes means of divergence time, standard error of means, 95% of highest posterior density intervals (HPD), effective sample sizes (ESS) 
and posterior probability of each clades. 



associated with Symbiodinium clades A, C or A + C combinations. In 
case of samples from the Indian Ocean Stylophora clades associated 
with Symbiodinium clades C, D or C + D. However, Stylophora 
samples from Yemen (clade 3) associated with Symbiodinium clade 
A and C similar to their clade 4 counterparts. This is also clear from 
the results of phylogenetic analysis in which the samples form Yemen 
(clade 3) form their own separate sub-clade. Differential 



Symbiodinium association could be due the presence of corals in 
different regional environments, which can significantly influence 
the ecology and evolution of this relationship 6 . 

The COI phylogeny of S. pistillata is corroborated by corallite 
structure (Fig. 2B), although extensive morphological variation still 
exists 7 . Generally, corallites are round in shape and regularly distrib- 
uted along the branches. However, corallite arrangement along the 
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branches varies among the 4 clades (Fig. 2 for detailed description of 
the skeletal properties in 4 clades). 

Divergence time between main clades was inferred using relaxed 
molecular clock based on COI sequences obtained from this study. 
Molecular clock calculation was calibrated using the earliest fossil 
record of Stylophora (S. octophylla, also see 'Divergence time' section 
in the Material and Methods for explanation of fossil record), found 
in Oman and dated around 65-70 Mya 8 . Mean divergence time of 
clades 1 and 2 from clades 3 and 4 is estimated to be > 51.46 mya 
(node A), followed by the divergence of clades 3 and 4 from each 
other (~34.55 mya, node C), the split of clades 1 and 2 from 
Seriatopora (—41.65 mya, node B), and finally the divergence of 
clades 1 and 2 from each other (—28.77 mya, node D) (Fig. 3). 
These changes correspond to significant changes in marine connec- 
tivity during the middle Eocene (54-38 mya) and Oligocene (38- 
23 mya) 5,9,10 in which the circumtropical Tethys seaway began to 
fragment, leading to the isolation of Europe from the Indian 
Ocean, and the relocation of the centre of diversity for shallow mar- 
ine species (including corals), from the west Tethys to the Middle 
East/East- African region 11,12 . Fragmentation of the Tethys seaway 
may therefore have driven the isolation of S. pistillata and the emer- 
gence of these different clades 12,13 . 

Discussion 

These results highlight a number of issues in the evolution and tax- 
onomy of this very common coral species. First, our results support 
the view that, for some marine taxa, the highest biodiversity is 
recorded in the Indian Ocean 10 , particularly the seas around 
Arabia, rather than the west Pacific. While three distinct clades over- 
lap along the Red Sea, Arabian Peninsula, African coast, Madagascar, 
and the oceanic islands in the west-central Indian Ocean, the entire 
Pacific Ocean is dominated by a single clade. Second, the genus 
Stylophora appears to be paraphyletic, as Seriatopora hystrixl calien- 
drum are embedded at the base of the PWA and CMSA clades. This 
supports similar findings for many other conventional taxonomic 
groups of corals, and further indicates that the present classification 
of corals, even at the family and genus level, needs major revision 1416 . 
Morphological convergence, phenotypic plasticity, recent speciation, 
and hybridization are likely causes of such patterns/phylogenies. 
Third, we show that the Stylophora "pistillata" complex consists of 
at least four clades. Corroborating this with a detailed morphological 
and microstructural characterization will be an important step 
towards formally resolving the taxonomy of this common coral spe- 
cies. In the Gulf of Aden, a recent genetic study 4 discriminated 
between two Stylophora pistillata clades, with a third clade including 
specimens from the western Pacific 4 . Moreover, morphologic and 
morphometric data showed the clades to be two distinct species 4 , 
Stylophora pistillata and Stylophora madagascarensis Veron, 2000. 
In the latter, strong genetic divergence was noted between the Indian 
and western Pacific populations, despite no difference in skeletal 
morphology. We suspect that two of the four clades evident in the 
present study could correspond to the above-mentioned species, with 
the other two clades comprising undescribed Stylophora species. 
These results demonstrate that the widely used binomen Stylo- 
phora pistillata actually encompasses several identities. 

Four distinct clades of S. pistillata highlight uncertainty in both the 
taxonomy 3,4 and systematics of corals 14,15 . This inhibits potential 
action at the species level 16,17 to conserve corals and coral reefs, such 
as assisted migration 18 which depends on a definitive understanding 
of biogeographic ranges in order to avoid unintended species intro- 
ductions, as well as policy action to identify species which are crit- 
ically endangered 19 . Multiple identities in S. pistillata cast doubt on 
prior research based on a single species 7 , and challenge the paradigm 
of coral as phenotypically plastic organisms with limited genetic 
variation 20,21 . 



Methods 

Collection of coral samples. Stylophora pistillata samples was collected from 34 
locations in reefs of Kochi and Okinawa (Japan), Tioman Island (Malaysia, West 
Pacific), Indonesia, Lord Howe-Australia, Western Australia, Lizard Island- 
Australia, Taiwan, New Caledonia, Chagos archipelago, Mauritius, Reunion, 
Madagascar (North and South), South Africa, Zanzibar, Kenya, Dijbouti, Yemen, 
Oman, Eilat (Israel, Red Sea), Saudi Arabia (Arabian Gulf, Gulf of Aqaba) and Iran. A 
small fragment of coral was clipped from each colony, placed in a labeled polyethylene 
bag, and preserved in 70% (v/w) ethanol. 

Molecular analyses. DNA extraction was conducted using the protocols described 
previously 22 . Genetic analyses of coral host were using two mitochondrial DNA 
region (Mitochondrial cytochrome oxidase I and tRNAW) and two nuclear DNA 
regions (ITS and Calmodulin intron 3). Sytnbiodinum clades of Sytlophora hosts were 
identified using two ribosomal DNA regions. Finally, the molecular divergence time 
of Sytlophora was estimated using relaxed molecular clock method. 

Mitochondrial cytochrome oxidase I (COI). A fragment of COI was initially 
amplified from S. pistillata using the primers, LCO1490 and HC02198, of 23 . PCR 
procedure was following 24 ; the products were directly sequenced from both ends 
using the same primers. All the sequencing was carried out by Mission Biotech Co., 
Ltd, Taipei, Taiwan. The sequences obtained were deposited in the Dryad Repository: 
http://dx.doi.org/dryad.n2fb2. 

ITS rDNA. Targeted segments containing the ITS1-5.8S-ITS2 region were amplified 
using the "anthozoan-universal" primer pairs, IS: 5'-GGTACCCTTTGTACAC- 
ACCGCCCGTCG CT-3' and 2SS: 5 ' -GCTTTGGGCGGCAGTCCCAAGCA- 
ACCCGACTC-3 ' , as described in 25 . PCR was performed using the following thermal 
cycle: 1 cycle of 95°C for 4 min; 4 cycles of 94°C for 30 s, 50°C for 1 min, and 72°C 
for 2 min; and 30 cycles of 94 C for 30 s, 55°C for 1 min, and 72"C for 2 min. The 
amplification reaction used 50 — 200 ng of template DNA and BRL Taq polymerase 
in a 50 ul reaction volume, using the buffer supplied with the enzyme, under 
conditions recommended by the manufacturer. The PCR products were 
electrophoresed in a 1% agarose (FMC Bioproduct, Rockland, ME, USA) gel in 1 X 
TAE buffer to assess the yield. PCR products were cloned using the pGEM-T system 
(Promega, Madison, WI, USA) according to the manufacturer's recommendations. 
Nucleotide sequences were determined for complementary strands of at least 
3 — 5 clones from each sample. The sequences obtained were deposited in the Dryad 
Repository: http://dx.doi.org/dryad.n2fb2. 

Calmodulin intron 3. Primers for calmodulin intron 3 (CAD3) were developed in the 
laboratory of C. A. Chen. A calmodulin-encoded cDNA clone was isolated from the 
egg-derived cDNA library of Acropora muricata 26 . This cDNA clone was 
characterized and compared to calmodulin genes available in GenBank. One primer 
was designed specifically for S. pistillata, CADStyl-F (5 ' -ACAAATGAAGTT- 
GGTGCTGATGGTAGGAGT-3')(Chiu et aL unpublished data) and the universal 
primer, CADCORAL-R (5'CTCTGGGAAGTC AATAGTGCCGTTTCC-3 ' ). PCR, 
cloning, and DNA sequencing were conducted following the conditions described for 
ITS. The length of CAD3 from different locations was compared. The sequences 
obtained were deposited in the Dryad Repository: http://dx.doi.org/dryad.n2fb2. 

Phylogenetic analyses. COI, ITS1 and CAD3 DNA sequences were aligned using 
ClustalW algorithm implemented in MEGA 5.0 27 . Alignments were then manually 
adjusted to increase the overall similarity. COI sequences (546 bp) were aligned and 
translated to amino acid codes using MEGA 5.0 across all samples. A deletion of ca. 
400 bp in one group of CAD3 sequence that was present from 177 bp to 593 bp was 
adjusted to final 380 bp after alignment. The ITS1 sequences with —1200 bp length 
were used in the analysis. 

Phylogenetic analyses of COI, ITS and CAD3 based on the Bayesian analysis were 
performed using MrBayes 28 and Maximum-Likelihood (ML) using PhyML 3.0 29 
(http://www.atgc-montpellier.fr/phyml/). For Bayesian analysis, the optimal 
molecular evolution model was determined using the Akaike Information Criterion 
(AIC) 30 performed in MrModeltest 2.3 31 . The most suitable models selected by the 
AIC for the analysis were HKY+I, GTR+I + G, and HKY for COI, ITS, and CAD3 
respectively. Markov Chain Monte Carlo (1,000,000 cycles) parameters used to run 
MrBayes were; 1 cold chain and 3 heated chains with 2000 for Burnin. Model selection 
for ML was run using the program jModeLTest2 32 within time-reversible (GTR) 
model under AIC. Models selected were TrN with rate variation among sites TrN + G 
(nCat— 4, gamma shape parameter— 0.5960) for COI, GTR+G (gamma shaper 
parameter— 0.5310) for ITS, and HKY (transition/transvertion ratio — 1.9428) for 
CAD3. The levels of robustness were assessed by 1000 bootstrap replicates for COI 
and CAD3, 500 for ITS for both Bayesian and ML analyses. Only COI sequences were 
analyzed using Neighbor-joining (NJ) (performed using MEGA 5.0), with 1000 
bootstrap and Kimura p-distance parameter was selected as the substitution model. 

tRNAW. A section of mtDNA, from NAD5 gene to partial COI gene 33 was amplified 
with primer sets FNAD5.2deg (5'-GCCYAGRGGTGTTGTTCAAT-3') and 
RCOBdeg (5'-CGCAGAAAGCTCBARTCGTA-3') and followed 34 . One of the 
genes, tRNAW, was found to have duplicated in some of the regions. PCR products 
were sent for direct sequencing then compared for duplication. To clarify different 
patterns for tRNAW in 4 different Stylophora clades, representative DNA samples 
from 4 clades were amplified using the primer set FNAD5.2deg and RCOBdeg and 
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were subjected to restriction enzyme digestion using Rsal (Thermo Scientific, UK). 
The reaction mixture was according to the manufacturers instructions. Amplified 
PCR products were incubated overnight at 37 C with RSAI enzyme. The digested 
samples were then run on 3% agarose gel {1% normal agarose + 2% low melting 
agarose) for 3 hours at 50V. The restriction pattern on the gel was photographed 
using gel documentation system (Vilber Lourmat, France) and visually analysed. 

Molecular phylotyping of the Symbiodinium clades. Molecular phylotyping of the 
Symbiodinium clades was conducted using 28S rDNA region and ITS2-DGGE. The 
protocol for 28S rDNA was modified from methods described by 35 . The variable 
domains, Dl and D2, of the 28S rDNA of Symbiodinium were PCR-amplified using a 
host-excluding primer set of Dl (5 ' -CCCGCTGAATTT AAGCATATAAGTA- 
AGCGG-3') and D2 ( 5 ' - GTTAG ACTCCTTGGTCCGTGTTTC AAG A-3'), and 
then digested with the restriction enzyme, Rsa I, to produce restricted fragment length 
polymorphism (RFLP) patterns, followed the protocol of 36 . The All enzymes were 
purchased from MBI (Fermantas). ITS2 region were amplified using primer set 
TTSintfor2' (5'-AAT TGC AGA ACT CCG TG-3') and TTS2 clamp' (5'-CGC CCG 
CCG CGC CCC GCG CCC GTC CCG CCG CCC CCG CCC GGG ATC CAT ATG 
CTT AAG TTC AGC GGG T-3 ' ) from 36 and using touch-down PCR 31 . PCR products 
of ITS2 were electrophoresed using 45-80% denaturing gradient gels for 16h on CBS 
Scientific system (Del Mar, CA, USA). Gels were stained with 1 X SYBR Gold (Life 
Technologies, Invitrogen, U.S.A.) for 20 min, and were photographed using gel 
documentation system (Vilber Lourmat, France). Method used for assigning the 
ITS2-DGGE fingerprint was followed as per 36 . Prominent bands of each fingerprint 
were sent for direct sequencing then match with the sequences from Genbank. 

p-distance analysis. The p- distances (number of base substitution per site ± standard 
error estimate) based on COI sequences within and between the four clades were 
determined using MEGA5 27 under Kimura 2-parameter model 37 . Interspecific 
pairwise distances in Stylophora and Seriatopora were compared with distances in 
other Scleractinia genera. COI sequences obtained in the present study were used for 
Stylophora and Seriatopora and COI sequences for other genera were obtained from 
the NCBI GenBank database. 

Divergence time. Divergence time of each clade was calculated (based on 
phylogenetic COI analysis) using Beast 38 , which allow relaxed molecular clock among 
different lineages. Yule birthrate process was chosen as prior, the distribution of 
divergence on each node was set as normal distribution with 5% of standard error. 
HKY model (Hasegawa, Kishino, Yano 85 model) with proportion of invariance 
(HKY + I) was selected as the most appropriate evolutionary model in the evaluation 
of likelihood ratio tests for molecular clock estimate 39 . 10 million generations were 
performed and saved every 1000 generations to calculate their phylogenetic 
relationship. First quarter of 10000 topologies were discarded as burnin while the 
remaining were saved to calculate the posterior probabilities. Fossils of Stylophora can 
be traced back to late Cretaceous in both the Caribbean and the Indo-Pacific region 
before it disappeared in the Caribbean during the early Miocene 40 . It appears that 
Stylophora genus was distributed worldwide prior to Cenozoic before becoming 
extinct in the Caribbean. To give a reliable reference point for molecular dating, 
wherein present day Stylophora occurs only in the Indo-Pacific, we used the 
information from the fossil record found in Oman (Stylophora octophylla), which can 
be traced back to its first appearance in Santonian and in Oman during Maastrichtian 
which is around 65.5-70 Mya 8 . 

Morphology analysis. Microstructure of coral skeleton from Taiwan, Chagos, 
Madagascar, Red Sea and Yemen was compared using photographs obtained through 
scanning electron microscopy 4 . 
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